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A primordial magnetic field in the early universe will cause Faraday rotation of the linear po- 
larization of the cosmic microwave background generated via Compton scattering at the surface of 
last scattering. This rotation induces a non-zero parity-odd (B-mode) polarization component. The 
Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data puts an upper limit on the magni- 
tude of the B-polarization power spectrum; assuming that the B-polarization signal is totally due to 
the Faraday rotation effect, the upper limits on the comoving amplitude of a primordial stochastic 
magnetic field range from 6 x 10~* to 2 x 10~® G on a comoving length scale of 1 Mpc, depending 
on the power spectrum of the magnetic field. 

PACS numbers: 



The cosmic microwave background temperature and polarization anisotropy measurements presented in the 5-year 
data of the Wilkinson Microwave Anisotropy Probe (WMAP) provide the most comprehensive characterization to 
date of the perturbations present in the early universe at the time of photon-baryon decoupling p[] . The microwave 
background polarization field is conventionally decomposed into a sum of parity-even ("E-type") and parity-odd ("B- 
type" ) parts; since B-type polarization is absent in a cosmological scenario involving only scalar density perturbations, 
it serves as a powerful test of other cosmological physics, including primordial tensor and vector metric perturbations 
[2]. Another source of B-type polarization is Faraday rotation of the orientation of linear polarization due to the 
presence of a cosmological magnetic field: if a purely E-type polarization pattern is everywhere rotated by 7r/4, it 
becomes a purely B-type polarization. Smaller rotation angles generate a component of B-type polarization from an 
initial pure E-type polarization field. 

Faraday rotation of polarized radiation is in general a powerful tool to estimate magnetic field amplitudes and 
spectra (see [3] for a recent review). In the cosmological context, Kosowsky and Loeb |1| proposed microwave 
background Faraday rotation to probe a homogeneous primordial magnetic field; this issue has been re-addressed by 
several authors for the homogeneous field [51|B], as well as for a stochastic field [TJ [51 HJ HD] • For a homogeneous field, 
the same polarization rotation effect additionally induces non-zero parity-odd cross correlations between microwave 
temperature and B-polarization anisotropies and E-polarization and B-polarizations anisotropics and correlations 
between different multipoles [12] which are absent for standard cosmological scenarios. Primordial magnetic fields are 
interesting because they could serve as seed fields for the observed magnetic fields in galaxies and galaxy clusters [TS] . 

Similar effects (microwave background birefringence and parity-odd anisotropy spectra) may also have non-magnetic 
explanations. In particular, microwave background polarization plane rotation occurs in models with violation of parity 
symmetry and/or Lorentz invariance in the early universe |14| . Parity-odd polarization anisotropy spectra appear 
due to gravitational lensing [15] or non-zero primordial magnetic helicity |16j . 

Here we use WMAP 5-year B-polarization power spectrum upper limits [T7| to place limits from Faraday rotation 
on a stochastic primordial magnetic field. We assume that the cosmological magnetic field was generated during or 
prior to the early radiation-dominated epoch. The high conductivity of the primordial plasma results in a "frozen-in" 
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magnetic field, fixing the temporal dependence to be the simple scaling B(x,r7) = B(x)/a^ where a is the scale factor 
and J] conformal time. Throughout this paper, B represents the comoving value of the magnetic field. We also 
normalize the scale factor a by setting qq — 1 today, and we employ cosmological units with h = I = c and gaussian 
units for electromagnetic quantities. 



II. CMB FARADAY ROTATION EFFECT 



Electromagnetic waves propagating in a magnetized medium have the plane of their linear polarization rotated (see, 
e.g., [18j). A linearly polarized wave can be expressed as a superposition of left and right circularly polarized waves. 
The magnetic field induces a phase velocity difference between the two circular polarization, resulting in rotation of 
the polarization plane. The rotation angle a for a plane wave with comoving frequency vq propagating in the direction 
n satisfies 

^ = U \2 2 " ■ B(x) (1) 

where f — Xenef^TO. is the differential optical depth, rie and Xe are the comoving electron number density and ionization 
fraction, ut is the Thomson cross section, and q is the magnitude of the electron charge. Faraday rotation of the 
microwave background is a subtle problem, because the polarization is generated and rotated simultaneously in the 
region of the last scattering surface. However, as long as the total rotation is small compared to 7r/2, the total rotated 
polarization angle can be expressed simply as an average of the rotated polarization angle from each infinitesimal 
piece of path length through the surface of last scattering, neglecting depolarization effects [8j . 

When considering a stochastic magnetic field, we make the simplifying approximation that any magnetic field 
component with a wavelength shorter than the thickness of the surface of last scattering is neglected: for these 
components, the rotation of polarization generated at different optical depths will tend to cancel, leaving little net 
rotation. Then we can treat the magnetic field as constant throughout the rotation region, so that the total rotation, 
which is the sum over the rotations of each infinitesimal piece of generated polarization, can be expressed as the total 
rotation incurred by the polarization generated at some particular effective optical depth (for details see [S].) 

A Gaussian random magnetic field is described by the two-point correlation function in wavenumber space as 

(S;(k)i?,(k')) - {2nf6^'\\. - k')P,,PB{k), (2) 

where Psik) is the magnetic field power spectrum, vanishing for all wavenumber larger than the damping scale 
ko, Pij = 5ij — kikj is the plane projector, and ki = ki/k. We use the Fourier transform convention i3j(k) = 
/ (Pxe^^'^Bjiyi). (We neglect nonzero magnetic helicity since it does not affect the polarization plane rotation [71 
[HI [IS].) The power spectrum Psik) is related to the total energy density Eb of the magnetic field through Eb = 
(PkPB{k), and it is given by simple power law PB{k) = ABk"'^. Smoothing the field on a comoving scale 
X > Xd = ^Tr/ko hy convolving with a Gaussian smoothing kernel, the smoothed magnetic field amplitude Bx is [S] 

i.^-^3, -^."(W2 + 3/2) ^ ^^^^ 
We assume that the magnetic field cut-off scale is determined by the Alfvcn wave damping scale [2U] . 



Mpc^V yiO-^Gj VMpc 

which will always be a much smaller scale than the Silk damping scale (thickness of the last scattering surface) for 
standard cosmological models. 

For a given comoving radiation frequency vq, consider the rotation angle Q;(n) of the microwave background linear 
polarization as a function of sky direction n. For a stochastic magnetic field, it is obvious that the rotation angle 
averaged over the sky is zero. Expanding the rotation angle in spherical harmonics leads to the definition of the 
angular power spectrum C" via 

9/ I 1 

(a(n)a(n')) = ^ ^CfFKn • n'). (5) 



where Pi{x) are Legendre polynomials. The rotation angle power spectrum due to a stochastic magnetic field is [8] 

91(1 + 1) Bl /^\"«+' r 

(4^)3(7^1.4 r(,,^/2 + 3/2) Uoj Jo 
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FIG. 1: The solid line shows the B-polarization power spectrum due to Faraday rotation from the WMAP best-fit cosmology, 
plus a stochastic magnetic field with amplitude Bx — 0.2 fiG, A = 1 Mpc, i/ = 30 GHz, and power law index ub = —2.9. The 
data points show binned WMAP 5-year B-polarization data with I < 800. Note that this magnetic field amplitude and power 
spectrum is ruled out by the data in the region between I = 300 and I = 500. Points with I < 150 do not contribute significantly 
to the constraint. 



where Xd = korjo, are spherical Bessel functions. This rotation multipole expression contains a sharp short- 

wavelength cutoff ku in the magnetic field; in reality, the effective cutoff will be smoothly spread over a range of 
scales. Given this rotation angle power spectrum, we need to compute the B-polarization power spectrum induced by 
the rotation field from the primordial E-polarization. This is given by Eq. (44) of Ref. [S], 

^BB ^ N?Y. ^^\+'^^f + NfK{lMMfCfi^Cr, {C\%^,f (7) 

where Cf^^i^Q are Clebsch-Gordan coefficients (with techniques for numerical evaluation in Appendix B 
of Ref. [H]), the normalization factor Ni = {2{l - 2)!/(/ + 2)!)i/2 and the function K{l,li,l2) = 
-i (L^ j^Ll + Ll- 2L1L2 - 2LiL + 2Li - 2L2 - 2L) with L = l{l + 1), Li = h{li + 1), and L2 = hih + !)• 
The cross-correlations between temperature and B-polarization and between E and B-polarization are zero for a 
stochastic magnetic field with zero helicity. The modifications to the existing E-polarization power spectrum and the 
cross-correlation between E-polarization and temperature are negligible for small rotation angles, and we ignore this 
second-order effect in our analysis. 



III. RESULTS AND DISCUSSION 



To compare the data with the theory prediction, we calculate C; for B\ = 1 nG, A = 1 Mpc, ub between -2.9 
and -1 in steps of 0.005, and / values from 2 to 1000. Our range of spectral indices is based on plausible magnetic 
field generation mechanisms [21 and spectral indices which are not excluded by other bounds; see Refs. [HJ [53]. 
The input Cf ^ values up to ^ = 4000 were calculated using CMBFAST for the best-fit WMAP cosmological 
model. The obtained set of is then further rescaled as i?^ and A"^+^ to obtain theoretical for different 

values of magnetic field amplitude and scale, respectively. WMAP data [5S] covers five frequency bands with central 
frequencies ranging from 23 GHz to 94 GHz. We use the measured B-polarization power spectrum limits for the three 
highest-frequency WMAP bands centered at 41, 61, and 94 GHz (Q, V, and W); our limits could be improved by 
including information from the lowest two frequency channels, but these are dominated by synchrotron foregrounds. 
We use only data for I > 32, for which we can treat measurement errors for different I values as uncorrelated. Most of 
the constraining power of the data comes from multipoles between I = 200 and I = 500, as displayed in Fig. 1, where 
the theoretical signal begins to be comparable to the errors on the data points; neglecting the data with I < 32 only 
marginally affects the limits obtained here. 

The theoretical Cf ^ power spectra are compared with the WMAP data using statistics. We conservatively 
include an extra 40% error on theoretical C^^ values to compensate for using the nominal frequency center of each 
WMAP band rather than the detailed frequency response of each band. Theoretical uncertainties on calculated values 
of Cf^ axe neglected, since these are much smaller than the experimental uncertainties on WMAP data. As the data 
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is consistent with the hypothesis of C^^ — (see Fig. [ij), we set limits on values of B\ as a function of ub and A. 
The 68% and 95% confidence limit bands on Bx as a function of are given in Fig. [2] We also display in Fig. |3]the 
95% confidence upper limits on Bx for fixed values of ub ranging from -2.9 to -1.0 as functions of A. 

These results presume that B-polarization signal is totally due to the Faraday rotation effect. Other possible 
source of B-polarization, such as inflationary gravitational waves gravitational lensing [T5 or non-zero primordial 
magnetic helicity [16] , will strengthen the magnetic field limits derived here. 

A primordial magnetic field on a scale A = 100 Mpc must have an amplitude Bx of less than 10^^ G for any 
power spectrum between ub — —2.9 and ub = —1. These limits are weaker that those obtained recently in Refs. |10j . 
Improvements on this limit using the Faraday rotation signal will be challenging, requiring substantially more sensitive 
polarization measurements at frequencies of 50 GHz or below; most current or planned high-sensitivity polarization 
experiments use bolometer detection technology, which tend to lose sensitivity in this frequency range. Upcoming 
measurements from the LFI instrument on Planck will be of roughly comparable polarization sensitivity to the WMAP 
measurements at similar frequencies. One notable effort is the QUIET experiment [29 , which uses coherent detector 
technology at frequencies of 40 and 90 GHz. It is eventually anticipated to measure the BB polarization power 
spectrum with an improvement in sensitivity over WMAP by a factor of 10^. The lowest frequency of 40 GHz, 
compared to 23 GHz for WMAP, makes the Faraday polarization power spectrum smaller by a factor of 7.3 for the 
same magnetic field; such an experiment could thus place limits on the amplitude of a primordial field which are a 
factor of approximately (10'*/7.3)^/^ = 37 more stringent than those here. A detailed analysis for this experiment 
[30j projects limits on a homogeneous primordial field below 10^^*^ G; however, primordial field limits will always 
be stronger than the stochastic field considered here, because Faraday rotation of a primordial field also generates 
nonzero TB and EB polarization cross-correlations which are larger than the corresponding BB power spectrum. 
Polarized foreground emission is also a potentially difficult systematic limit to these measurements. Other upcoming 
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polarization experiments aimed at detecting the BB polarization from inflation are likely to improve on the Faraday 
rotation limits here, though the Faraday rotation power spectrum amplitude decreases by a factor of 230 between 
WMAP's 23 GHz channel and a likely 90 GHz lowest-frequency bolometer channel. 

The cosmological magnetic field itself generates microwave background anisotropics, particularly via the tensor and 
vector perturbations it induces (see Refs. |23J for overviews of magnetized cosmological perturbations). Here we 
have neglected these perturbations, limiting ourself to considering the Faraday rotation effect alone, which provides 
a distinctive signature of magnetic fields nearly independent of other properties of the universe. We use the direct 
observational data without any priors. Our magnetic field limits are weaker than those arising from CMB temperature 
maps [ini HZ] or current bounds on a homogeneous magnetic field due to correlations between different I modes [25] 
and CMB temperature non-gaussianity |31j . The strongest future limits on magnetic fields in the early universe will 
likely come from limits on the temperature and polarization fiuctuations generated by magnetic field-induced vector 
perturbations [26J; we will address these current limits from WMAP data elsewhere. 
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